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An Effective Solution to the Nonlinear, Nonstationary 
Navier-Stokes Equations for Two Dimensions 

By Ralph E. Gabrielsen, Ames Research Center and U.S. Army Air Mobility R&D 
Laboratory, Moffett Field, California, USA 

1 . Introduction 

Classical fluid flow problems of aerodynamics stem from the nonlinear, 
nonstationary Navier-Stokes equations* When solving these problems for cases 
in which simplification is unacceptable, the practice has been to proceed in 
one of two ways: (1) the problem is considered within the framework ol the 

Reynolds equations, or (2) the Navier-Stokes equations are dealt with numeri- 
cally. There are, however, significant disadvantages to both the Reynolds 
equations approach and to existing numerical solutions. Relative to (1), 
recent experimental evidence indicates essential information about the flow 
dynamics is lost in the process of generating the Reynolds equations 
[1, p. 325], and relative to (2), existing numerical solutions of the Navier- 
Stokes equations completely lack explicit error estimates. 

Whereas the development of more effective techniques for solving such 
fluid flow problems by means of the Reynolds equations is highly unlikely, 
future progress seems to lie in the direction of improved methods of handling 
such problems numerically. 

Towards this end, a numerical solution with explicit error estimates is 
presented for the Navier-Stokes equations in two dimensions with suitable 
boundary conditions. 
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The solution, based on Newton's method as generalized by Kantorovich and 
Akilov [2], [3], Is the limit of a sequence of functions that are solutions to 
a sequence, of related linear partial differential equations. 

Background material, preliminary results, principal development of 
the solution, and conclusions are presented in the sections that follow. 


2. Background 

The nonlinear, nonstationary Navier-Stokes problem is mathematically 
described by 


- VAAl|l - + fly “ f2jj I 

H.OS) - bj , AiK3S) - b2 , - ♦(S) , ) 

where S Is the Interior of a compact two-dimensional region, and 3S its 
boundary. 

It will be assumed that » where 

C,j.^(S) = {f(x,y,t)|N times continuously differentiable in (x,y) 6 S,T times 
continuously differentiable in t,0 < t < » , inclusive of spatial derivatives 
up through order [N/2]} . Note that [ ] denotes greatest integer part. 

C^^CS) will be considered under the norm 


11^1 


m 

T L?J n 
T 1-0 n»0 m-o 


.nti, 

3 

3tSx"ay"-"' 


0<t<®» 



max 
x.y^s 
m-o 0<t<» 


3”i|; 

3x”3y"-"“ 


Let 


P(l|)) = At|)^ - vAAi|/ - 'I'yAll/jj + 'I'x^'l'y + fZjt " ^ly • 
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Equation (1) can be expressed in the equivalent form: 

P(i|i) - 0 , iK3S) - bi . AiKaS) - bj , 

P will be considered a mapping from Ci**(S) into Co*^(S). 

3. Preliminary Results 

Lemma 1. The Frechet derivative of P exists, and is given by 
P'(«o) - A - vAA - «|.0yA - AH»o^ ^ + D^q^A ^ + AtJ»0y • 

Proof* The Frechet derivative is defined as follows: 

||P(i|; + a) - P(4») - La|l^ o 

if n — n “ 0 for some linear operator L, then 

f|a|ho 

- L . Hence it can be directly shown that for <|»q € Cj**(S), P'('I'q) 
exists, and 

F'<n) - a ^ - *0/ ^ - *0, iy * ^ • 

(See ref. [4].) 

Lemma 2. P''(ifr) exists for all i|» € Cj‘*(S), and 

P"('l<A)4»e - d. A6 + 0 A«ji - 0 A(|) - (fi A0„ . 

^y X y ^x X y x y 

In addition, 

l|P"('l'o)ll i 1 ’ 

Proof. By definition, if there exists a bilinear operator B such that 

lln + -''‘V 

then P"(()»jj) exists and P"(i|»q) - B. Hence it follows directly that 

P"(.Pq)* 0 - <^yA0^ + 0yA^^ - 0^A^y - (>^A0y 

and that ||P"||^1[4]. 
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In order to show that the "Newton-Kantorovich Algorithm" 
0 -P0)-„) +P'(V('('„rt-l - '»'m> 


ij.j^(9S) - bj , - b2 I (3) 

'f'mit-o " ’ m - 0, 1, 2, , I 

converges to a solution of Eq. (2), it is sufficient to show [2], [3] that 

l|p(4»o)ll l|p’(>Po)’MP l|p"('<'o)ll • 

Since for a given II^(4 ^q/II Is readily calculated and 

||P**(4>o)|| ^ If it remains to calculate an upper estimate for I 1*^ * I I • 

4 • Main Development 

Since, for a given ^ Cj^(S), the sequence of arising from Eq. (3) 

converges to the solution of Eq. (1), providing 

\\f(%)\\ l|p'(’i-orMP l|p"(i'o)M ‘ 

it remains to determine if and when P'(i{|q)~^ exists, and estimate 

To this end, let m = 0 in Eq. (3). Therefore, 

0 - P (iPq) + P • (<!»(,) (i|<i - iI/q) in CqO ; 

ij»jj(9S) ■ bj , <J»i(9S) “ bj , A\J»q( 9S) « bj , Ai|/j(9S) » b2 , 

’*'o|t-0 “ '''l|t-0 “ • 

Let ^ therefore, -P(<|»q) • P'C'J'q)^* ♦ exists, then 

^ » -P’-i(t|»o)P(t|»o) . 

Hence, we seek a solution i|i of the problem 

0 ■ POI/q) +P’(i|'q)^ , Mas) » 0 , A|(|(9S) » 0 , and - 0 . (4) 

Let 

P’(4'o> - A ^ - uAA + A , 
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where A is a linear operator mapping into Cq^(S). 

“ ^'J'Ox^y * ^'*'0y^x 

Equation (4) can therefore be expressed by 

di|)^ - vAAi - -P(<l»o) - A(i|(jj)iji , 5(3S) • 0 , A 4 I 0 S) “ 0 , ^jt»0 " ° 

Let iji be defined as follows: 


-If II 


[A(i|>n)iJi](x",y",T) 


+ P(i|»o)(x".y".T) } <7 (x.y.x'ty'.x",y",t - T)dx" dy" dx dx' dy' , 


where 


C'’(x,y,x' ,y* ,x",y",t) = G(x,y,x* ,y' ) h'^(x' , y* ,x",y",t) . 
Lemma 3. If i £<]• (^) exists, then i satisfies Eq. (4). 
Proof. Since 

C'^(x,y,x' ,y' ,x",y",t) « G(x,y,x' ,y' )h'’(x' ,y* ,x”,y",t) 
where G is the Green's Function for La Place's eouation, and h'’ the 
kernel function for the heat equation 

(j)j. - vA(() = 0 , <|'(3S) = 0 , '*’|t=0 “ ® * 

it readily follows that 

i|»OS) = 0 , At|»(3S) * 0 , and “ 0 • 


Also, 

A^(x 


I. 7. 

J. 


A'|>j.(x,y,t) - - 


[A(ij»n)|(x",y",T) + P(«J/.)(x",y",T)]H'’(x,y,x",y",t - x)dx" dy" dx 


[A(t|»Q)|(x",y",x) +P(i|»jj)]H^''(x,y,x",y",t- x)dx" dy" dx 


[A(<>-Q)|(x",y"t)+P(.|«o)(x",y",t)]H'’(x,y,x",y",0)dx" dy" . 

AA|(x,y,x)--| I [A4 )i(x",y",x)+P(4.o)(x",y",x)]AH'’(x,y,x",y",t- .)dx« dy" dx 
•'O •'S 
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Therefore » 




[A(<|/q) +P(«(»Q)][H^'^(x,y,x'',y'',t - t) - vAH'^]dx" dy" dt 


-L 


[A()|»Q)i + P(i|/o)]H'’(x,y,x",y'',0)dx" dy 


■-A(l|)g)i - P('l'o) • 


Suf violent conditions under which 


m 

il> 


of Eq. (6) exists will now be 


developed. Define the linear mapping B: 



q(x'',y",T)i3'’(x,y,x',y’,x",y",t- x)dx' dy' dt 


dx" dy" . 


Therefore, Eq. (6) c.n be expressed in the form 

\p - BA\p - , 

(I - BA>;j - . 

Let Cg^(S) H the completion of Cg^(S) under the norm 

Ikll - max |t))| *f max |i|» I + max |i() I 
x,y€S x.yes * x,y€S ^ 

0<t<® 0<t<» Oit<» 


+ max |i(»^^| 

+ max |t(»^„| + max |i|/ 

x,y€S 

x,y€S x,y€S ^ 

0<t<*» 

0<t<® 0<t«» 

+ max |Ai(» 1 

+ max |Ai|; 1 • 

x,yes 

x,y€S y 

0<t<® 

G<t<® 


If (BA) is considered a mapping from Cg^(S) Into Cq^(S), and 

||BA|| < 1, then (I - BA)~^ exists, and 1(1 ■ (I - BA)“lW* (i(»g) . Therefore, 

In order to show that t|» exl~ 3, It Is sufficient to show that 

||B|| I |A| I < 1. To this end consider A as a mapping from Cq^ Into 

For ^ € C°, I |i()| L 0 ■ kl ; 

^0 x,y€S 
0<t«» 
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for ♦€ Ikllr 3 “ »ax kl + “** l*xl “* 


x,y€S x,y€S 

0<t«» 0<t<» 


*’ x,y€S y 
Ost<« 


+ max I + max |i|» | + max |i(i I 

x.y6S ^ x,y€S ** x.yes ^ 
0<t<« 0<t<* Ost<" 


Since A 


+ max I + max |Ai|» | . 
x,y€S * x,y€S ^ 

0<t<« Oit<» 

||A^|| by definition. 


sup 

lUIL -1 

Co^ 


A(i|»o)^ - ’*’o/*y ^'"Oy'^x ’ 


and 


||A^| L 0 • »ax |A^| « max |-<l»o - Ailfj, * + 'I'q + ^'I'Oy'^xl * 

0 x,y€S ^ X X 

0<t<» 

it follows that | |A| | < M, , where 

M - max (kovl * I’I'o * 

'"o x,y€s y “x X y 

0<t<“ 

Consider B as a mapping from into | |b| | - sup (||Bz||~3) 

° Ilz|lr0-1 0 


by definition. Therefore, since 


Bz| ~ 3 - max 


|Bz| + max |(Bz)jj|+max | (Bz) | + max U®*) 3 j 3 {l + ®®x | (Bz) 


xy ' 


B :s 


+ max UBz)yyl + max | (ABz)^| + max |(ABz)y| , 

max I I I H''(x',y',x",y",t-T)[G(x,y,x’,y') + |Gy| + Ig*! 
x,y€S •'S i'o •'S 


0<t<" 


+ IGxxI + |Gyy| + |Gxy|]dx" dy" dx dx' dy' 

+ max f 1 [|Hx''(x,y,x",y",t-t)| + |Hy''(x,y ,x",y",t - x) | ]dx" dy 

•'0 i^S 


'• dx 
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Let 


,y' t)(G 

5 


N. • max 
x,yf^S 
0<t<» 

+ lG^l)dx” dy” dT , 


I I + I Gy I + I Gjqj I + I Gyy | 


and 


M - max(|i|»o I, |if»Q I, |Ai|; |, |Ai|»q |) . 

Vq y ”x "x y 


Therefore, for M, N- < 1, i|» exists, and t|» ■ (I - BA)“^ BP(4 *q). Hence, from 

Vo 

Lemma 3, it follows directly that -P'(Vn) ^ • (I - BA)"^ B and that 






1 - IIBII l|A| 


- 1 - N«M, 

® Vo 


5. Conclusicn 

From the preceding sections it directly follows that if i|»q, the initial 
educated guess, is judiciously chosen, a sequence of functions 
i|»jj(n ■ 1, 2, . . .) can be directly constructed that converge to the exact 
solution of the nonlinear, nonstationary Navier-Stokes problem as expressed 
by Eq. (1). In addition, an explicit error estimate can be directly deter- 
mined for each approximate solution i|^^. 

Expllcl^y then, for 



with M, N- < 1, the sequence of ij^-'s determined by Eq. (3) exists and 
Vq * ” 

converges to the exact solution \li* of Eq. (1). 

Moreover, under these conditions, tue method yields the explicit error 
estimate 
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